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Abstract. In this paper, we present results of our calculations on the effects of next to nearest 
neighbor boson hopping (t') energy on Bose-Einstein condensation in cubic lattices. We consider 
both non-interacting and repulsively interacting bosons moving in the lowest Bloch band. The inter- 
acting bosons are studied making use of the Bogoliubov method. We find that the bose condensation 
temperature is enhanced with increasing t' for bosons in a simple cubic (sc) lattice and decreases for 
bosons in body-centered cubic (bcc) and face-centered cubic (fee) lattices. We also find that interac- 
tion induced depletion of the condensate is reduced for bosons in a sc lattice while is enhanced for 
bosons in bcc and fee lattices. 
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1. Introduction 

Studies of Bose-Einstein condensation in optical lattices and crystalline lattices is an active 
field of research both in atomic [1^1] and condensed matter physics [5-11]. In condensed 
matter physics, there have been extensive studies of bose condensation of bipolarons [5], 
exctions [6], exciton-polaritons [7,8], and magnons [9-11]. Studies of bosons in optical 
lattices may be said to have received a boost with the demonstration of bose condensed to 
Mott insulator transition [2] predicted in theoretical studies [1,12-15] of strongly interact- 
ing lattice bosons. In the presently available optical lattices, it has been shown [1] that it is 
sufficient to include nearest-neighbor (NN) hopping of bosons in the kinetic energy part of 
the Hamiltonian of the system. Nevertheless, considering the fast pace of developments in 
this field, it may be useful to investigate the effects of the next to nearest-neighbor (NNN) 
hopping on bose condensation in optical and crystalline lattices. Recently, we presented a 
study of the lattice symmetry effects on bose condensation in cubic lattices [16]. In that 
work, we had confined to NN hopping of lattice bosons. In this paper, we extend this 
work including NNN boson hopping. The bosons are considered to be of spin-zero and 
charged (see also the note in Ref. [17]). We would like to emphasize that we are not exclu- 
sively considering bosons in optical lattices. Our calculations should be considered in the 
enlarged context including bose condensation in crystalline lattices. In the next section, 
we describe the models and methods used in our calculations along with a discussion of 
results. The conclusions are given Sec. III. 
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2. Bose condensation in cubic lattices with NNN hopping 

Non- interacting bosons: Consider bosons moving in cubic lattices. The energy eigen- 
functions of a single boson moving in a periodic optical or crystalline lattice potential are 
Bloch waves [18] and energy eigen- values form bands. The Hamiltonian of non-interacting 
bosons in an energy band is: 

= 5>(fc) - M]4cfc , (1) 
fc 

where e(k) is the one-boson energy band structure, k is the boson quasi-momentum, /i 
is the chemical potential, and ct is a boson creation operator. Within a tight-binding ap- 
proximation scheme [19], including the NN and the NNN Wannier functions overlaps, the 
i-band structures we consider for cubic lattices are: 
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where the lattice constant has been set to unity. Here t is NN boson hopping energy and t' 
is NNN boson hopping energy in the lattice. 

The condensation temperature Tb for bosons in these bands can be calculated from the 
boson number equation 

11 ~ AT AT AT X/X/X] c{k x ,k y ,k^)- li > (5) 

1 *x ly y 1 ' z , , , „ k^T 1 

" ka, k y kz e B — 1 

where N s = N x N y N z is the total number of lattice sites, ks is the Boltzmann constant, 
T is the temperature, and n is number of bosons per site. We have numerically solved the 
bosons number equation (Eq. 5) to obtain the bose condensation temperature and ground 
state occupancy. The results of these calculations for various lattices considered are shown 
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Figure 1. The bose condensation temperature vs NNN hopping t' for noninteracting 
bosons in various cubic lattices. The curves are shown for: sc (solid line), bcc (dotted 
line), and fee (dash-dot line). These results are for n =0.25 (top panel), n=0.4 (middle 
panel), n=1.5 (bottom panel). In this and other figures W is the half-band-width. 
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Figure 2. The variation of condensate fraction with temperature (T) for bosons in a sc 
lattice with NNN hopping t' : t' = (solid line), t' = t/10 (dotted line). Here n =0.25 
(top panel), n=0.4 (middle panel), n=1.5 (bottom panel). 

in Fig. 1. We find that t' increases the bose condensation temperature of bosons in a sc 
lattice. For bosons in bcc lattice the Tb decreases with increasing t'. For bosons in a fee lat- 
tice also increasing t' more or less leave Tb unaltered. These trends can be approximately 
understood in the low boson density limit. In this limit, states with significant thermal pop- 
ulation is close to the bottom of the energy bands. Now, for a given small boson density, 
the ratio between bose condensation temperature and half-band-width is proportional to 
l/(m*W), where m* is the boson effective mass. We find that, in the low density limit, 
Tb/W (where Tb is the bose condensation temperature and W is the half -band-width) 
goes as: (t + At') /3f for sc, (t + t')/{2t + 3t') for bcc, and 1 /3 for fee. On plotting, one 
can easily see that the Tb /W increases with t' in the sc case, decreases slightly for bcc 
case, and remains constant for the fee case. These trends are consistent with our numerical 
results. 

The growth of the condensate fraction and the number dependence of Tb for bosons in 
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Figure 3. The bose condensation temperature vs n for bosons in a sc lattice with NNN 
hopping t' : t' =0 (solid line), t' = t/20 (dash-dot line), t'=t/10 (dotted line). 



a sc lattice shown in Fig. 2 and Fig. 3 are similar to that found for the case of t' 
Similar results are obtained (not shown) for bosons in bcc and fee lattices. 
Interacting bosons: The Hamiltonian of interacting bosons is 



H 



fj]c\c k 



U 
2Nl 



[16]. 



(6) 



where U is a constant boson-boson repulsive interaction energy. To treat the effect of 
interactions we make use of Bologliubov approach [20] to the interacting bosons system. 
In this theory, it is assumed that the ground state of interacting bosons system is a bose 
condensate. Since the condensation occurs into the lowest single particle state (for which 
k = in our cases), one gets < ejeo > ~ < cqc\ >. This allows one to treat the operators 
and Co as complex numbers and one gets < cj > = < cq > = ^Nq. Here Nq = noN s 
in which no is the boson occupancy per lattice site in the k = state. The second order 



interaction terms are obtained from the substitution: c - 
approach, the boson number equation is obtained to be (for details see Ref. 16): 
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where £(fc) = e(k) — cq, £o is me energy of the lowest single particle state, and E{k) = 
v / e(k) + 2Un ^(k). The primes on the summation signs indicates that the sums excludes 
the k — state in to which the bosons condense. The Bologliubov method used would 
valid so long as the interaction energy is smaller than the kinetic energy of the bosons. 
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This approximately translates to U < 2W. The effect of increasing interaction (U) is to 
lead to an increase in the effective mass of the bosons which eventually gets localized for 
large interaction strengths . But, this happens for integer filling. Our results (Figs. 4-6) 
are for n = 0.25, 0.4, and 1.5 which is not close to integer filling. The interaction induced 
enhancement of the boson effective mass will not be significant in this case since there are 
sufficient number of un-occupied sites in the lattice so that the bosons can move around 
without paying a penalty for multiple boson site occupancies. The condensate fraction for 
bosons in various cubic lattices are shown in Fig. 4-6. 

For bosons in a sc lattice, we find that interaction induced depletion of the condensate 
is reduced with increasing t' as shown in Fig. 4. For bosons in a bcc lattice (Fig. 5), 
increasing t' is found increase the interaction induced depletion. In the case of bosons in 
an fee lattice, the effects of increasing t' does not have much effect on condensate fraction 
as shown in Fig. 6. The bose condensation temperature is unaffected by the interaction in 
the Bogoliubov method. 

3. Conclusions 

In this paper, we investigated the effects of NNN hopping of non-interacting and interact- 
ing bosons in cubic lattices on bose condensation temperature and ground state occupancy. 
We find that the bose condensation temperature is enhanced with increasing t' for bosons 
in a simple cubic (sc) lattice and decreases for bosons in body-centered cubic (bcc) and 
face-centered cubic (fee) lattices. We also find that interaction induced depletion of the 
condensate is reduced for bosons in a sc lattice while it is enhanced for bosons in bcc 
and fee lattices. These results would be relevant to bosons in condensed matter systems 
in which NNN boson hopping is not negligible. The results could also be applicable to 
bosons in optical lattices. For instance, it was recently shown that hard-core lattice bosons 
moving in an optical lattice and interacting with phonon modes of polar molecules trapped 
in the lattice develops significant NNN hopping amplitudes [21]. There is a hope that sev- 
eral models of strongly correlated quantum many particle systems can be simulated in a 
controlled manner in optical lattice systems [3,4]. As mentioned earlier, there is also a 
renewed effort in investigations of bose condensation in crystalline lattices. Further, one 
of the routes to superconductivity is through the condensation of charged bosons (bipo- 
larons, for example). Furthermore, higher temperature superconductivity may be possible 
in correlated electron systems by the condensation of charged bosons generated within 
the electrons system through strong correlation effects [22,23]. Though we are unable to 
find a concrete example at this juncture, it is not inconceivable that the energy spectra of 
some of the possible emergent boson modes in strongly correlated lattice electrons systems 
disperses away in momentum space with significant contributions from NNN hopping am- 
plitudes. Then, our results may have some relevance to these bosons as well. 
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Figure 4. The variation of Condensate fraction (for T=0) with U/W of weakly inter- 
acting bosons in a sc lattice with NNN hopping t'\ t' = (solid line), t' = t/20 (dash-dot 
line), t' = t/10 (dotted). Here n =0.25 (top panel), n=0.4 (middle panel), n=1.5 (bottom 
panel). 
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Figure 5. The variation of Condensate fraction (for T=0 and n=0.4) with U/W of 
weakly interacting bosons in a bcc lattice with NNN hopping t' : t' = (solid line), t' 
= t/20 (dash-dot line), t' = t/10 (dotted line). Here n =0.25 (top panel), n=0.4 (middle 
panel), n=1.5 (bottom panel). 
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Figure 6. The variation of condensate fraction (for T=0 and n=0.4) with UAV of 
weakly interacting bosons in a fee lattice with NNN hopping t' : t' = (solid line), 
t' = t/20 (dash-dot line), t' = t/10 (dotted line). Here n =0.25 (top panel), n=0.4 (middle 
panel), n=1.5 (bottom panel). 
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